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1. Introduction 

The study of triangular-lattice spin-| Heisenberg model (THM) has been a central 
problem in quantum many body physics ever since Anderson made the proposal that 
its ground state properties could be described by a Resonant Valence Bond picture 
[TJ [2] . The development of several numerical techniques [21 IH El El El El El HOI ED] was 
crucial to elucidate that the zero point quantum fluctuations, in this particular system, 
are not enough to destroy the classical 120° Neel order, lending support to a simple 
semiclassical description at low temperatures, such as that provided by linear spin wave 
theory (LSWT) [J2]. 

Nonetheless, early high temperature expansion (HTE) studies [6] performed down 
to low temperatures showed no evidence of a renormalized classical behaviour as 
predicted by the non linear sigma model (NLoM) [131 EH]. For instance, around 
T = 0.25 J, the correlation length calculated by HTE is only £ ~ 1.5 lattice constants, 
in contrast to the value £ ~ 12 predicted by NLaM [SI [TO]. Consistent with these 
values, the entropy calculated by HTE was one order of magnitude larger than that 
of the NLaM. These early results were interpreted as a probable crossover between 
renormalized classical and quantum critical regimes |15| . 

Unexpectedly, series expansion (SE) studies [T0l[T6] performed at zero temperature 
also showed a strong downward renormalization of the high energy part of the spectrum 
with respect to LSWT, along with the appearance of roton-like minima at the midpoints 
of the edges of the hexagonal Brillouin zone (BZ). Originally, the presence of such roton- 
like excitations were proposed to be related to possible fermionic spinon excitations 
which in turn would lead to the anomalous low temperature properties of the THM 
[TO] . However, subsequent works showed that non trivial 1/S corrections, arising 
from the non collinearity of the 120° Neel order, accurately recovers the T = series 
expansion results [J71 [181 [19]. This gave support to an interacting magnon picture 
for the spectrum, although it was found that the magnon-quasiparticles are not well 
defined for a significant part of the Brillouin zone [T81IT9]. Nevertheless, by assuming a 
bosonic character for the SE dispersion relation, the high values of entropy found at low 
temperature was attributed to the thermal excitation of rotons, even at temperatures 
below the roton gap [TO] . 

Here we explore the low temperature properties of THM from an alternative 
viewpoint: A bosonic spinon perspective, based on the Schwinger boson formalism 
[201 |2T] . One advantage of this point of view is that it preserves the rotational invariance 
of the THM at finite temperatures, in agreement with the Mermin- Wagner theorem [22J ; 
while at zero temperature it recovers the 120° Neel ordered state as a condensation of 
the Schwinger bosons [231 121] • In fact, it has been shown in the literature [251 126] that 
the Schwinger boson theory describes very well the ground state properties of the THM, 
finding good agreement with exact diagonalization [27] and other available numerical 
techniques (For a complete survey of the available results on the THM see table III 
of [10]). More recently, we have shown [28] that the main features of the magnetic 
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excitation spectrum found using series expansions can be reproduced correctly using a 
proper mean field scheme of decoupling. In addition, by computing the density-density 
and spin-spin dynamical structure factors we were able to identify unphysical magnetic 
excitations that can be traced back to the relaxation of the local constraint of the 
Schwinger bosons at the mean field level |28j. 

In practice, the relaxation of the local constraint seems to be not crucial for 
the correct description of certain ground state static properties [28J. But as soon 
as temperature increases, the system starts to explore in increasingly large amounts 
an unphysical phase space, leading to an incorrect estimation of the thermodynamic 
quantities. Instead of introducing ad hoc factors to compensate for the above problem 
[20| |2T| . we use our ability to distinguish between the physical and the spurious 
excitations to properly extend the Schwinger boson mean field (SBMF) to finite 
temperatures by considering only the physical low energy excitations. We rename the 
latter the reconstructed Schwinger boson mean field (RSBMF). The values of entropy 
and uniform susceptibility thus calculated interpolate quite well within the temperature 
range — 0.3 J, on opposite ends of which LSWT plus 1/5 corrections and HTE become 
reliable, respectively. Our results support the idea of [TO] that the high values of entropy 
found with HTE are due to the excitation of rotons which, within the context of our 
theory, can be identified with collinear short range AF fluctuations above the underlying 
120° Neel correlations. 

2. Rotational invariant Schwinger boson mean field theory 

In using Schwinger boson representation for the spin operators [20] 



s? = \$lfa + b\M sf = jfil.k, - 5^), Sf = \{h^ - blki), (i) 

the local constraint on the boson number J7 ct b ia bi a = 2S must be imposed to fulfil the 
spin-S" algebra. The spin-spin interaction of the triangular-lattice Heisenberg model can 
then be written in terms of singlet bond operators [2"B"| |2"§] as 

Sj-Sj = : B^Bij : —A^Aij, (2) 

where Ajj = | J2 a ab ia bj^ and Z?J- = \ J^ CT b io .bj a are singlet operators, invariant under 
SU(2) transformations of the spinor (6^,64) and :: means normal order. The 
identities B\-Bij = 2{Si + Sj) 2 and A|^4y = 2(Sj — Sj) 2 reveal the ferromagnetic and 
antiferromagnetic character of each term of ([TBI , respectively. Hence, the possible 
coexistence of both kinds of correlations renders this scheme of calculation ideal to 
investigate frustrated quantum antiferromagnets [251 ESI I2H1 1221 GEH ED] • 

The mean field decoupling of f lT3|) is performed in such way that A^ = (Aj,) = (Aj) 
and Bij=(Bij) = (Bj-) . Then, after introducing the local constraint on average through 
a Lagrange multiplier A, the diagonalized mean field Hamiltonian becomes 



H M f = -Egs + ^ w k d^Akt + & ] _ H a-ki 



(3) 
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Here 

E_ = 
2 



£ S s = ^ 5>»< + AiV(S + i) (4) 

k 

is the ground state energy and 

wkt = <M = u k = [( 7k B + A) 2 - ( 7 £) 2 ]^ (5) 

is the spinon dispersion relation with geometrical factors 7^ = ^J^2§B$ cos k. 5 and 
7k = |</ ^5 As sink. 5, an d the sums run over all the vectors 5 connecting the first 
neighbors of the triangular lattice with N sites. The mean field parameters have been 
chosen real and satisfy the relations Bs = B_s and A$ = — A_s [28]. The mean field free 
energy is given by 

F = £ gs + T^m(l-e-^), (6) 

ko- 

whose minimization with respect to the mean field parameters, As, Bs and A leads to 
the following self-consistent equations: 

As =^E — (l + 2n k )sink.<5, 



2N ^ u k 
k k 



B* 



1 ^M±^) ( n_ 2nk)cosk . 5) (7 ) 



2N ^ co k 
k 

2N ^ u k 



2 2N ^ w k v ; 

k K 

with n^ = (e^ 0Jk — 1) _1 the bosonic occupation number. As previously stated, the 

advantage of the SBMF is to be able to study finite temperature rotationally invariant 

phase of the triangular AF, dictated by the Mermin- Wagner theorem [20J EI]. In 

particular, the numerical self consistent solutions of (J7]) correspond to the renormalized 

classical regime with an exponential decay of the spin correlation functions characterized 

by the magnetic wave vector Q = (^, 0) of the 120° Neel structure [32]. This manifests 

in the gapped spinon dispersion u-^ with minimum at k = ±y. The finite temperature 

gap prevents the infrared divergences from appearing in the theory. At low temperatures 

the gap is exponentially small, becoming u,q ~ 1/N in the T — y limit. Formally, one 

2 
should first perform iV — y 00 and then take the limit T — y to recover the SU(2) broken 

symmetry state. In analogy with the Bose condensation phenomena, this procedure 

requires one to treat separately the singular modes of (J7J) and transform the sums 

into integrals. This gives a new set of equations with the quantum corrected local 

magnetization m as a new self consistent parameter, while A is adjusted so as to get a 

gapless spinon dispersion in each iteration. Then, antiferromagnetism is interpreted as 

a quantum fluid where the condensate of the up/down bosons at ±y and the normal 

fluid of bosons correspond to the spiraling magnetization m and the zero point quantum 

fluctuations, respectively [21]. Alternatively, we have worked with large finite systems. 

This procedure is simpler because the same set (j7]) are used for the T^O and T = 
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cases. Although the mean field solutions keep their rotational invariant character, the 
local magnetization m can be obtained by relating it to the static structure factor 
SQ*) — SR elk ' R (s s l'S'o - 5'R|gs) evaluated at k = Q and the singular modes as 



(7l + A) 2 N 

S(Q) = —m 



2 



2N wi v " / 2 



where the last line of (JTj) has been used. For the triangular S — | AF the local 
magnetization of the 120° Neel ground state gives m = 0.275 (We have checked that 
this procedure is equivalent to solve the new set of equations above mentioned with self 
consistent parameters As, Bs and m and the sums transformed into integrals). This 
value should be compared with the quantum Monte Carlo result [7], mQMc — 0.205(1). 
The underestimate of the zero point quantum fluctuations can be attributed to both 
the mean field approximation and the relaxation of the local constraint which violates 
the physical Hilbert space. In the next section, we will show that the latter leads to 
some unphysical magnetic excitations in the spectrum, requiring a careful analysis of 
the low lying energy excitations to properly compute the low temperature properties of 
the THM. 

3. SBMF low lying magnetic excitations 

To analyze the magnetic excitations of the THM within the SBMF approximation we 
compute the T = spin-spin dynamic structure factor, 

S zz (k,oo) = Y, \(gs\Sl\n)\ 2 5(uj - (c - E 9 )), (8) 

n 

where \n) are mean field excited states, and S£ is the Fourier transforms of Sf. We have 
recently shown that (jHJ) takes the simple form [28] 

S ZZ (k, Ul) = ^ I -»k+q^q ~ ^q^q+k| 2 5(u; - (w_ q + Cd k+q )) , (9) 

q 

where u^ = [|(1 + — — )]a and fk = isgn(7^)[|(— 1 + — — )]a are the coefficients of the 
Bogoliubov transformation that diagonalizes ^mf- It is clear from (|9]) that the spin-spin 
dynamical structure factor consists of two spinon excitations which gives a continuum 
of spin-1 excitations. Furthermore, since we are working with finite systems whose 
mean field solution does not break the SU(2) symmetry, S xx = S yy = S zz . Nevertheless, 
as the system size increases, and the 120° Neel correlations are developed, there is a 
distinction among the two-spinon processes that can be observed in figure [T^a), where 
an intensity plot of ()9]) in energy- momentum space is shown. Notice that, in order to 
show the contribution at all energies, the intensity is plotted in a logarithmic scale. The 
main signal of S zz (k, uS) (intense red curves) corresponds to the microscopic processes 
of destroying one spinon & ± q ct of the condensate and creating another one b q in 
the normal fluid and viceversa; while the blue region corresponds to the creation of 
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Figure 1. Dynamical structure factors along the path shown in the inset of (b). 
(a) Magnetic structure factor within the Schwinger boson mean field © (Intensity 
curves); LSWT results (solid green line); series expansion results [TO] (blue dots). 
The reconstructed dispersion relation ZJk is shown in dashed line, (b) Density-density 
dynamical structure factor (fTTj) . The intensity scale is logarithmic as in (a) 



two spinons in the normal fluid only [28 J. The energy cost of the former two spin-| 
spinon excitation is just w Wrr Q, since in the thermodynamic limit uj^q — > 0. Notice 

K + 2 + 2 

that at low energies these shifted spinon based dispersions reproduce quite well the 
LSWT results (solid green line) which describe the semiclassical long range transverse 
distortions of the 120° Neel order. On the other hand, at higher energies, where the 
LSWT semiclassical description is no longer valid, the spectral weight between both 
shifted spinonic bands is redistributed in such a way that if one reconstructs a new 
dispersion cJk from those pieces of spinon dispersions with dominant spectral weight 
(dashed line), the main features of the series expansion results pU] (blue points) are 
recovered. There is, however, an additional remnant weak signal which we attribute 
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to the relaxation of the local constraint of the mean field approximation. To quantify 
the relative spectral weight between oJk and the weak band, note that in the best case 
-at k = Q- the contribution of cJk represents around 95% of the total weight while the 
weak band only 2%. On the other hand, in the worst case - middle point of OC - the 
contributions are 50% and 20%, respectively. To describe the physical Hilbert space of 
the spin operators the local constraint of Schwinger bosons must be satisfied exactly, 
Sf = y( y + 1). So, no fluctuations of the number of boson per site should be allowed. 
However, as the constraint is taken into account on average there appear unphysical spin 
fluctuations in S(k, u) coming from such density fluctuations which can be investigated 
through the density- density dynamical structure factor defined as, 

N(k,u) = J2 |(gs|n k |n)| 2 % - (e n - E gs )), (10) 

n 

where n^ is the Fourier transform of rii = ^ a b\j)i a . After a little of algebra it results 
in 

N(k,Uj) = 22 l U k+q' U q + ^q^q+k| 2 5(w - (w_ q + W k +q))- (H) 

q 

In figure [H(b), (TTTj) is plotted in energy-momentum space. It is observed that the 
dominant spectral weight of the density- density dynamical structure factor iV(k, to) 
coincides with the weak signal of S zz (k, u) (figure Ufa)). This notable resemblance 
led us to identify the remnant signal of S zz (k, u) with the unphysical spin fluctuations 
originating from the density fluctuations of the Schwinger bosons [2H]. Therefore, we 
expect that after projecting the mean field ground state |gs) into the constrained Hilbert 
space of 25 bosons per site such unphysical excitations will disappear. 

So far the reconstructed dispersion corresponds to a spin-1 excitation made of two 
spin- 1 free spinons. Then, the relevant question is whether the spinons are bound 
or not once one goes beyond mean field theory. This issue has been addressed for 
frustrated AF within the context of effective gauge field theories [33j US]- I n particular, 
for a commensurate spinon condensed phase, it is expected that the fluctuations of 
the emergent gauge fields confine the spinon excitations, giving rise to the spin-1 
magnon of the 120° Neel order. While this kind of calculation is out of the scope 
of our present work, we can get some indication of the above physical picture by 
a simple calculation. Let us split the original Hamiltonian as H = Hmf + V, with 
the interaction term given by V — H — Hmf- The effect of V on a two free spinon 
state 1 2s) = |qa;pcr) = a^at |gs) can be estimated, to first order in perturbation 
theory, by computing the energy of creating two spinons above the ground state as 
-E , 2s = (2s|if |2s)— {gs\H |gs). Then, the interaction between the two spinons is obtained 
as v int = E 2s - E^ F where E^ F = (2s\H MF \2s) - (gs\H M F\gs)= ^V + u; pa and E 2s is 
rescaled as E 2s = \E 2s in order to compensate the difference, (gs\H |gs)=|(gs|iJMF|gs), 
resulting between our mean field decoupling and the application of Wick's theorem 
which corresponds to a fully self consistent Hartree-Fock-Bogoliubov decoupling 



(Note that this | has nothing to do with the ad hoc factor of Arovas and Auerbach 
). The interaction thus calculated give 
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1 



: , y L 7q-p(MqMp + l^ql>p| 2 + ^V^UpVpj + (12) 

+ 27q+p(«qM 2 + Mp|Wq| 2 - 2w q W q W p t> p J + 9 J , 

by the present mean field where 7k = |^ 5 J#cosk ■ 5. For a physical excitation 
1 2s) = | ± y°"i k "F 2 a )' i nv °l v i n g the creation of one spinon in the condensate and 
another one in the normal fluid, it is easy to check numerically the attractive character 
of (TT2"j) . In particular, when the total two spinon momentum is k =B the energy 
binding is t>i nt ~ — 0.16,7, while the energy cost to create the two free spinons above 
the ground state is cuq +w r _q ~ 0.7 J. On the other hand, when the two spinons are 

2 a 2 

created in the condensate the interaction becomes |t>i nt | ~ 0(Nm 2 ), meaning an infinite 
attraction for the spinons that build up the magnons at the Goldstone modes k = 0, ±Q. 
Even if this instability is an artifact of the first order correction, we believe that this 
simple calculation lends support to the physical picture of tightly bound spinons in 
the neighborhood of the Goldstone modes while at higher energies they remain weakly 
bound. Of course, this statement should be addressed by a more rigorous calculation 

ESI. 



4. Low temperature properties 

As pointed out in the introduction, the SBMF ground state properties like energy, 
magnetic wave vector, magnetization and spin stiffness are not significantly affected by 
the relaxation of the local constraint [251 123 EQ]. However, as soon as temperature 
increases the system starts to explore in increasingly large amounts an unphysical phase 
space due to the fluctuations of the density of Schwinger bosons. Consequently, the 
SBMF becomes inadequate to describe the low temperature properties. In this section 
we will show that it is possible to modify the SBMF in order to get reliable results for 
the low temperature regime of the THM. First, we assume that at the temperatures 
considered, the spinons are sufficiently bound in such way that the relevant physical 
excitations can be envisaged as spin-1 excitations with the reconstructed dispersion 
relation cJk defined above. On the other hand, in order to mimic the projection to 
the physical phase space -not considered in the mean field- we discard the unphysical 
excitations. We have renamed this calculational scheme as reconstructed Schwinger 
boson mean field (RSBMF). We will mainly focus on the importance of the roton-like 
excitations to explain the high values of entropy found at intermediate temperatures 
(T ~ 0.3 J) [HI OH]- Although the RSBMF solutions correspond to the renormalized 
classical regime, it is worth stressing that we do not expect to fully recover the 
expected behaviour near the zero temperature transition. In fact, it is well known 
that the correlation length at the mean field level has an extra exponent in temperature 
dependence of the preexponential factor when compared with the predictions of effective 
field theories based on the NLoM or confined spinons near the zero temperature 
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Figure 2. Temperature dependence of the reconstructed dispersion relation, as defined 
in the text, along the path of the inset of figure [1] The horizontal dashed lines 
correspond to the energy intervals used in figure 01 



transition J35J|36]- We have found non trivial solutions, As, Bs ^ 0, for the self consistent 
equations ([7]) up to T ~ 0.42J. For T > 0.42J the finite temperature phase of the 
triangular AF becomes a perfect paramagnet with no intersite correlations. Such a 
phase has no analog in the interacting spin problem, so it is an artifact of the mean field 
or large-N approximation. In fact, the effect of finite N corrections to this unwanted 
transition has been investigated in the literature [37]. In any case, since the available 
results of HTE are reliable down to T ~ 0.3 J we will concentrate on the temperature 
range — 0.42J in order to interpolate between other results. 

In figure [2] is shown the temperature dependence of the reconstructed dispersion 
oJk- As temperature increases, a narrowing of the bandwidth along with the opening of 
a gap and a flattening of the roton-like excitations is observed. The entropy per site 
corresponding to these low energy spin-1 excitations is 



S = — ^2 [ ( n k + 1) In (rik + 1) - n k In n k 

k 



(13) 



where the reconstructed dispersion Wk is plugged in the averaged occupation number 
l) -1 . In contrast, when the two free spinon species are considered (SBMF) 



n k 



/3u) k 



the resulting entropy per site is 

<S = — ^2 t (^ka + 1) In (n ka + 1) - n kCT In n ka 



(14) 



kcr 



where the spinon excitations u; ko - are plugged in the occupation number n ko - and the 
sum over the two spinon species is taken into account. In figure [3] it is shown that the 
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Figure 3. Temperature dependence of the entropy computed within the RSBMF (|T3|) 
(black solid line); SBMF (03) (red dashed line); HTE [BJ (orange dotted-dashed line); 
SE [10] (blue dotted line) and NLctM [TO] (green double dotted-dashed line). The black 
thin solid line represents the SBMF result divided by 2.785 (see text in section [5]). 



RSBMF entropy (black solid line) interpolates quite well between the expected values 
at zero temperature and T ~ 0.3 J predicted by HTE (orange dotted-dashed line); while 
the high values of the SBMF entropy (red dashed line) are due to the inclusion of 
the spurious excitations. To complement these results, also shown in figure [3] is the 
entropy corresponding to the low energy spin-1 excitations with the dispersion relation 
found with SE at zero temperature (blue dotted line) and that of the NLoM [ID] . The 
agreement with SE is very good at low temperatures, although the RSBMF entropy is 
better aligned with HTE at higher temperatures. This difference can be attributed to 
the fact that the SE entropy has been computed with the T = dispersion relation while 
in the RSBMF entropy the reconstructed dispersion relation is temperature dependent 
(see figure [2]) . 

To discern between the contribution of the magnonic low energy excitation and the 
high energy roton-like excitations we have computed the entropy for different ranges 
of energy [10] . The blue dashed line of figure H] represents the contribution of the low 
energy modes while the red dotted-dashed line represents the low energy plus the roton 
modes. In agreement with [10J it can be observed that the increasing contribution of the 
roton modes starts from T ~ 0.1J. As a reference we have included the entropy value 
predicted by the NLcrM at T = 0.3 J (black cross) which agrees reasonably with the low 
energy contribution of the RSBMF. The apparently counter-intuitive decreasing of the 
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Figure 4. Temperature dependence of the RSBMF entropy computed for different 
range of energies. The cross indicates the entropy value predicted by the non linear 
sigma model (NLcrM) at T = 0.3 J. 



partial entropies at larger temperatures is due to the temperature dependence of ujy_. 
In the case of a temperature independent dispersion [10] it is expected an increasing 
contribution for all energy intervals because of the increasing occupation number average 
for each mode. In our case, instead, the effect of temperature is to push up an important 
amount of modes outside the corresponding energy intervals, lowering its contribution 
to the partial entropy for T > 0.3J. For instance, the minimum of the dispersion at 
T = 0.4J is greater than 0.65J, then the contribution of this energy interval to the 
entropy is zero. Similarly, within the energy interval — 0.9J, the decreasing of the 
entropy at T = 0.4J is due to the fact that the rotonic part of the dispersion flattens 
and crosses the top of the interval — 0.9J. 

Taking into account the above mentioned we conclude that at least within the 
range 0.1J — 0.3J the contribution of the rotonic modes becomes relevant to the total 
entropy. In particular, these results confirm the idea that the higher values of entropy 
found with HTE around T ~ 0.3J can be attributed to the contribution of the roton 
excitations. Recently, we have related the roton-like excitations with AF collinear 
fluctuations above the 120° Neel order [28]. In fact, it is easy to check that by including 
frustrating interactions to second neighbors ( J2) the roton excitations soften, becoming 
the precursor of the transition to a collinear state at J2/J ~ 0.16 |39j. Then, within the 
context of the RSBMF, the high values of entropy at low temperatures can be related to 
the contribution of AF collinear fluctuations signalled by the flattening of the dispersion 
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relation around the midpoint of the edges of the triangular Brillouin zone. This effect is 
missed in the NLaM, where only the effective low energy modes are taken into account. 
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Figure 5. Temperature dependence of the uniform susceptibility. RSBMF (TT5j) f Solid 
black line) ; SBMF (flo) (red dashed line) ; HTE of [6] (orange dotted dashed line) and 
SBMF computed within the one singlet operator scheme as in [32j [40] (blue dotted 
line). The arrow indicates the averaged T = result of the linear spin wave theory 
plus 1/S corrections of [41] . Inset: temperature derivative of the uniform susceptibility 
versus temperature. 



The uniform susceptibility calculated within the RSBMF gives 

S(0) i r _. 



Xu 



T 



(n k + 1) 



(15) 



where S(0) is the finite temperature static structure factor S(q) = J2r < ^o-Sr > 
exp(iq.R) evaluated at q = 0; while in SBMF it results 



Xu = jj^2n k (n k + l). 



(16) 



In figure |5] is shown the uniform susceptibility calculated within the RSBMF (black 
solid line) and the SBMF (red dashed line). Consistent with the previous entropy re- 
sults, the uniform susceptibility is overestimated by the SBMF. That is, the presence 
of the spurious magnetic excitations enhances the response of the system to a uniform 
magnetic field. On the other hand, the RSBMF results interpolates quite well between 
the expected zero temperature value (see arrow) and HTE results of [6] (orange dotted- 
dashed line). In particular, the extrapolated zero temperature value is Xu ~ 0.072 
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which should be compared with Xu ~ 0.084, corresponding to the averaged spin wave 
results plus 1/5 corrections [15]. At very low temperatures the RSBMF behaves as 
Xu ~ 0.072 + 0.6T while the behaviour expected for classical renormalized regime [15] is 
Xu ~ 0.084 + 0.07T. The one order of magnitude in the slope is related to the fact that 
the true correlation length increases faster than the mean field solutions as zero tem- 
perature is approached (see [151 EE]). For T > 0.15J the uniform susceptibility differs 
from HTE results in contrast to the excellent agreement obtained for the entropy. While 
entropy depends on the correct counting of states, through the reconstructed dispersion 
Uk, the uniform susceptibility is related to the spectral weight distribution at k = 
which at these temperatures, probably, is not correctly taken into account by the mean 
field approximation. For instance, we have assumed that the pair of spinons building 
up the spin-1 excitations at k = are sufficiently bound. Therefore, the required low- 
ering of the uniform susceptibility to recover the HTE results would be consistent with 
a weaker bound spinon regime for T > 0.15J. Another source of the above mentioned 
discrepancy can be attributed to the omission of the continuum contribution to the 
uniform susceptibility. Notice, however, that despite the rounded peak of the HTE is 
broader than the RSBMF they are both located approximately at the same temperature 
position, T ~ 0.35 J [6]. This can be noticed in the inset of figure [5] where the tempera- 
ture derivative of Xu versus temperature is plotted. For completeness, we also show the 
SBMF prediction within the one singlet operator scheme (blue dotted line), previously 
performed in [321 [40]. In this case the interpolation is not good because this scheme of 
calculation fails to reproduce the low energy spectrum of the THM [28J . 

Finally, at low temperatures (below T ~ 0.1J), the specific heat behaves 
quadratically as C v ~ a{T/J) 2 , with a = 7.3 and a = 5.2 for SBMF and RSBMF, 
respectively. As a reference, the latter value should be compared with a = 5.3(2), 
resulting from an HTE interpolation method [42], specially developed for the low 
temperature specific heat behaviour; while a ~ 3.4 is found by plugging in the spin 
wave velocities of the THM in the Debye construction. 

5. Discussion 

Here we discuss the validity of the RSBMF developed to compute the low temperature 
properties of the THM. Originally, the THM was investigated within SBMF based 
on the one singlet operator (A?) scheme. As this scheme gives a wrong sum rule, 
rE ka S aa (k,oj)duJ = %NS(S + 1), it is well known [201 E2] that an ad hoc factor of 
2/3 is needed to compensate the overcounting of the number of degrees of freedom. 
For the present two singlet scheme there is no need of ad hoc factors since already 
at the mean field level (SBMF) the sum rule is fulfilled. Furthermore, it has been 
shown that Gaussian corrections for the ground state energy and spin stiffness of the 
THM improve significantly the accuracy of the SBMF results when compared with exact 
diagonalization results on finite systems 
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Regarding the magnetic excitations, we have defined a reconstructed dispersion w k - 
guided by the dominant spectral weight of the spin-spin dynamic structure factor- which 
reproduces quite well the series expansion results. In addition, we found unphysical 
excitations related to the fluctuations of the density of Schwinger bosons due to 
the relaxation of the local constraint. Unfortunately, it is difficult to compute the 
Gaussian corrections for the dynamical structure factor to improve the local constraint. 
Consequently, to mimic the projection of the SBMF into the physical Hilbert space we 
have discarded the spurious excitations. This allowed us to define the reconstructed 
Schwinger boson mean field, based on bosonic spin-1 excitations with the reconstructed 
dispersion uJk- 

Given that the RSBMF thus defined is not rigorously justified we have followed an 
alternative route to confirm whether the counting of the degree of freedom is correctly 
taken into account, at least approximately, within the RSBMF. To do that we have 
concentrated on the entropy for the two-spin problem. In this case the problem can be 
solved exactly and in the high temperature limit the entropy is <S ex = 21n(2S* + 1). 
On the other hand, the SBMF prediction can be worked out analytically, giving 
S MF = 21n[(5 + l) 2S+2 /S 2S }. For S = 1/2, it is found that S MF = 2.785 S ex . Therefore, 
the factor 2.785 gives a faithful compensation of the SBMF entropy due to the wrong 
counting. Notably, if the SBMF entropy (red dashed line) of figure [3] is divided by the 
factor 2.785 (see black thin solid line of figure |3]) the RSBMF entropy (black solid line) 
of figure [3] is recovered, at least within the temperature range — 0.3 J of interest. This 
agreement gives further support to our RSBMF procedure. 

Another possible procedure is to implement the local constraint numerically by 
using a valence bond basis but sign problems, already at the variational level, appear 
in the Quantum Monte Carlo calculation [13]. More recently, a permanent algorithm 
was used to compute ground state energies on the kagome lattice within the context 
of projected Schwinger bosons. Although the method seems to circumvent the sign 
problem, only systems up to iV = 36 size has been studied [4"4"j . 

6. Concluding remarks 

We have investigated the low temperature properties of the triangular-lattice Heisenberg 
model with a bosonic spinon theory based on the Schwinger boson mean field theory. 
Using the two singlet operator scheme of calculation and by analyzing the spin-spin and 
density-density dynamical structure factors we can distinguish between the physical 
magnetic excitations, which reproduce the series expansion results, and the spurious 
excitations that result from the relaxation of the local constraint on the number of 
bosons per site. By assuming that i) the effect of projecting into the physical Hilbert 
space can be mimiced by getting rid of the latter excitations and ii) that an attractive 
residual interaction binds the spinons, we obtain a reconstructed Schwinger boson mean 
field with low energy spin-1 excitations with the reconstructed 'physical' dispersion 
relation. A comparison with reliable T = and high temperature (T > 0.3J) results 
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reveal that the RSBMF thus defined provides a very good interpolation for many 
thermodynamic properties like entropy, uniform susceptibility and specific heat over 
the temperature range < T < 0.3J, which is difficult to access by other methods. 
One of our main physical results is to confirm the idea [10] that the high values of 
entropy found with HTE are due to the excitation of roton excitations which, within the 
context of our theory, can be identified with the collinear short range AF fluctuations 
above the underlying 120° Neel correlations [28J. An alternative interpretation has 
been proposed in [15], using the physics of the XX Z Heisenberg model as a starting 
point. Here a low energy effective theory is constructed in terms of fermionized vortices 
where vortex-antivortex excitations on the honeycomb (dual) lattice are related to the 
roton excitations. Interestingly, the dependence of the roton excitations with spatial 
anisotropy resembles that with temperature found in the present work. Even if the 
description in terms of fermionized vortices recovers the roton excitations predicted by 
SE [TQl US], the present bosonic approach seems to be more appropriate to describe all 
the expected features of the isotropic AF Heisenberg model, that is, the 120° Neel order, 
the correct Goldstone mode structure and the roton excitations. In any case it would 
be interesting to investigate the validity range of each approach by performing a close 
comparison between them. 

The simplicity of the RSBMF along with the consistent description of several 
features of the THM like static ground state properties, energy spectrum and low 
temperature thermodynamic properties give strong support to the bosonic spinon 
hypothesis to interpret the physics of the triangular-lattice Heisenberg model. Of 
course, the present approximation should be refined by using many body diagrammatic 
techniques [35J or by performing 1/N corrections [30] . Work in this direction is in 
progress. 

Finally, we believe that a proper extension of the present theory to the anisotropic 
XX Z case would allow one to investigate the unusual magnetic features of the recently 
found [1HJ H3 HH] inorganic spin-1/2 triangular antiferromagnets, BasCoSb20g and 
Ba^CuSb20g. An experimental realization of the related XY model has been proposed 
recently in the area of ultra-cold atoms in optical lattice potentials [IH1 E0] • 
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